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Abstract 

We claim that the Higgsino-like and wino-like neutralinos can be good dark 
matter candidates if they are produced by the late time decay of Q-ball, which 
is generally formed in Affleck-Dine baryogenesis. The late time decays of the Q- 
balls into these LSP's and subsequent pair annihilations of the LSP's naturally 
lead to the desired mass density of dark matter. Furthermore, these dark matter 
can be much more easily detected by the dark-matter search experiments than the 
standard bino-like dark matter. 



Supersymmetry (SUSY) has been widely considered as an attractive framework for 
physics beyond the standard model. It explains the stability of the electroweak scale 
against quadratically divergent radiative corrections. Furthermore, particle contents of 
the minimal SUSY standard model (MSSM) lead to a beautiful unification of the three 
gauge coupling constants of the standard model. 

One of the remarkable features in the MSSM is the existence of an ideal dark matter 
candidate, that is, the lightest SUSY particle (LSP).[] The most extensively studied LSP 
as a dark matter candidate is the bino-like neutralino, since its thermal relic abundance 
naturally provides a desired amount of the present mass density of the dark matter. On 
the other hand, there have been much less interests in other candidates, such as the 
Higgsino-like and wino-like neutralino, since their thermal relic densities are generally 
too low to be a significant component of dark matter [I]. 

In this letter, we claim that the Higgsino-like and wino-like neutralinos are good 
dark matter candidates in spite of their large annihilation cross sections, if the origin 
of the observed baryon asymmetry lies in the Affleck-Dine (AD) baryogenesis 0, |3|. In 
AD mechanism, a linear combination of squark and/or slepton fields (0 field) has a large 
expectation value along a flat direction during an inflationary stage, and its subsequent 
coherent oscillation creates a large net baryon asymmetry. The coherent oscillation of 
the field is generally unstable with spatial perturbation and it fragments ||, [5|, ||] into 
non-topological solitons [0, Q-balls. It has been, in fact, shown in detailed numerical 
calculations || that almost all of the initial baryon asymmetry carried by the <fi field is 
absorbed into the Q-balls. 

The Q-ball has a long lifetime^ and its decay temperature is likely to be well below 
the freeze-out temperature of the LSP, which leads to the non-thermal production of 
the dark matter We show that in the case of the Higgsino-like or wino-like LSP 
the late-time Q-ball decay and subsequent pair annihilations of the LSP's naturally give 
rise to a desired dark-matter energy density. It is very encouraging that such candidates 
are much more easily detected than the standard bino-like neutralino. As we will see 
later, the detection rate is in fact more than ten times larger compared with the case of 



bino-like LSP flTJ. 



Let us first estimate the present energy density of the non-thermally produced LSP. 



1 We assume that the i?-parity is exact and hence the LSP is absolutely stable. 

2 We do not consider the gauge-mediated SUSY breaking Q, where the Q-ball is generally stable 0, 

!• 

3 In Ref. |6|, the authors only considered the case in which there is effectively no pair annihilation of 
the LSP after its production, and did not consider LSP's with large annihilation cross sections, such as 
Higgsino-like and wino-like neutralino. In general, their case leads to an overproduction of the LSP. See 
discussion below. (See also Refs. 
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Suppose that there is a non-thermal production of LSP at temperature T = T d , where T d 
is below the freeze-out temperature Tf of the LSP. (Tf is typically given by Tf ~ m x /20, 
where m x is the mass of the LSP.) The subsequent evolution of the number density n x 
of the LSP is described by the following Boltzmann equation: 



n x + 3Hn x 



— (av) n 



x ■ 



(1) 



where the overdot denotes a derivative with time, H is the Hubble parameter of the 
expanding universe, and (av) is the thermally averaged annihilation cross section of 
the LSP. Here, we have neglected the effect of the pair production of LSP's, which is 
suppressed by a Boltzmann factor exp(— m x /T) for T < m x . It is useful to rewrite the 
above equation in terms of temperature T and number density of LSP per comoving 
volume Y x = n x /s, where s = (27r 2 /45)g*T 3 is the entropy density and is the effective 
number of relativistic degrees of freedom. Here and hereafter, we assume that the energy 
density of the LSP is much smaller than that of the radiation p rad at T ~ T d , and no 
extra entropy production occurs after that. (We will justify this assumption later in the 
case of Q-ball decay.) Then, we obtainF] 
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where M p i = 2.4 x 10 GeV is the reduced Planck scale. This equation can be analytically 
solved by using approximations g*(T) ~ g*(T d ) ~ const and (av) (T) ~ const, which 
results in 

-i 



Y x {T d 
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45 



(av)M pl (T d -T) 
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Therefore, for sufficiently large initial abundance Y x (T d ), the final abundance Y x0 
T T d is given by 



Y, 



45 



(av) M pl T d 



for 



(4) 



Notice that the final abundance Y xQ is determined only by the temperature T d and the 
cross section (av), independently of the initial value Y x {T d ) as long as Y x {T d ) ^> Y x0 . 
From the above formula, we obtain the relic mass density of the LSP in the present 
universe: 
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4 Here, we have used g* p ~ <?* for T 
degrees of freedom for energy density. 
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1 MeV, where g* p (T) = (3Q/TT 2 )p rad (T)/T 4 is the effective 
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where h is the present Hubble parameter in units of 100 km sec _1 Mpc _1 and Q x = p x j p c . 
(p x and p c are the energy density of LSP and the critical energy density in the present 
universe, respectively.) Notice that the obtained abundance is much larger than the result 
for thermally produced LSP. In fact, it is enhanced by a factor of ~ (Tf/Td) compared 
with the case of standard thermal production with the s-wave dominant annihilations. 

Now let us discuss the LSP production by the Q-ball decay. First of all, the baryon 
number density ng is related to the number density of the Q-balls and the initial charge 
of each Q-ball Qf. 

Qi^q = fn B , (6) 

where / denotes the fraction of the total baryon asymmetry which is initially contained 
in the Q-balls. Notice that almost all the baryon asymmetry is initially stored in the 
Q-balls ||, namely, / ~ 1. 

The decay rate of a single Q-ball is given by |LH| : 



Q ~ ^~192vr 2 ' {> 



where u ~ m^, is the soft scalar mass of the <ft field, and A is the surface area of the 
Q-ball. (The radius of the Q-ball is given by R Q ~ y/2\K\- l / 2 m^ 1 , where \K\ ~ 0.01- 
0.1 0, H |15|.) Then, we obtain the lifetime of the Q-ball r d = Qi/Tq, or equivalently 
the decay temperature Td of the Q-ball: 

Actually, the formed Q-ball has a large charge and typically decays at T d ^ 1 GeV , in 
particular when the <fi field is lifted by a nonrenormalizable dimension-six operator in the 
superpotential with a cutoff scale ~ M p i. Hereafter, we will take T d ~ 10 MeV-1 GeV. 
From above equations, the production rate of the LSP per time per volume is given by 

N x T Q n Q = N x T Q ^-e(Q,-T Q t) 

= N x fn B « e -^±, (9) 

Td 

where N x is the number of LSP's produced per baryon number, which is at least 3. 
Thus, the evolution of the number density of the LSP is obtained by solving the following 
equation: 

n x + 3Hn x = NJ ( ?») sx - (av) n\ , (10) 



V S Jo T d 
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where we have normalized the baryon number density by the entropy density and sub- 
script denotes the present value. 

In Eq. fllOD, it is assumed that the LSP is uniformly distributed. Because the LSP's 
are produced from the Q-ball, which is a localized object, one might wonder if the pair 
annihilation rate of the LSP becomes much larger and its final number density becomes 
much smaller. However, we can see this is not the case as follows. First of all, it is found 
from Eq.(|]) that the number density of the LSP approaches its final value only after 
(Td — T)/Td ~ 0(1), which means it takes a time scale At ~ r^. (Notice that this is 
true for any local number density, as long as it is large enough.) By that time, LSP's 
have spread out by a random walk colliding with the background particles, and form a 
Gaussian distribution around the decaying Q-ball. The central region of this distribution 
has a radius f ~ ^/z/r^, where v~ x ~ G 2 F m x T^ ||. [Gp is the Fermi coupling constant.) 
Meanwhile, we can see that the number of Q-balls within this radius is much larger than 
one, roughly given by (4vr/3) f 3 n Q ~ 10 10 x (T d /1 GeV)- 6 (g i /10 20 )- 1 (m x /100 GeV)~ 3 / 2 . 
Hence, the assumption of the uniform distribution is justified. 

We should also note that the energy density of the Q-balls pQ is much smaller than 
that of the radiation p ra d for T ^ T d : 

p Q m^nQQi 
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where we have used the fact that the energy of the Q-ball per charge is roughly given by 
m^. Therefore, no significant entropy production takes place during the Q-ball decay. 

We have numerically solved the Boltzmann equation (|T0"D for Higgsino-like and wino- 



like neutralino and calculated the relic abundance f2 x , where we have taken N x = 3, 
/ = 1, (riB/s)o = 0.7 x 10~ 10 , and h = 0.7. In our calculation, we included final 
states; W+W~, ZZ, tt, h°A° } H°A°, Zh°, ZH°, and W ± H^^ (We took the cross 
sections from Ref. fl6f .) The results are shown in Figs. [l|-|3] in the m x -T^ plane. Figs. [I] 
and H correspond to the Higgsino-like LSP, where we took Mi = (3/2)/i and M 2 = 
3/i, while Fig. || corresponds to the wino-like LSP, where we took Mi = (3/2)p and 
M 2 = (l/2)/x. (Mi and M 2 are the soft gaugino masses for the SU(2)l and U(l)y 
gauge groups, respectively, and p denotes the SUSY contribution to the Higgs-boson 
(Higgsino) masses.) As for the other parameters, we used tan/3 = 5, m^o = 300 GeV 

5 Here, we included only s-wave annihilation cross sections, which is a reasonable approximation for 
Higgsino-like and wino-like neutralino and for T <^ m x . 
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and all trilinear scalar couplings a = 0, in all figures. For sfermions we assumed universal 
soft masses my. = m . We used m = 1 TeV in Figs.[l] and |3|, and m = 330 GeV in 
Fig.0. We used the g,{T) given in Ref. |7| for 100 MeV<T< 1 GeV, adopting the 



QCD phase transition near 150 MeV. 

It is found from these figures that both of the the mass densities of the Higgsino-like 
and wino-like LSP in fact fall in the desired region Q x ~ 0.1-1 in a wide range of LSP 
mass, for Q-ball decay temperatures ~ 10 MeV-1 GeV.[] (We have also confirmed 
that these results are well reproduced by the analytic calculation given in Eq. (|5|).) It is 
remarkable that the Higgsino-like and the wino-like LSP's can be excellent dark matter 
candidates even in the relatively small mass region, where the thermal production would 
give rise to too small relic abundance. (See discussion below Eq.@.) As we will see 
later, these regions are also advantageous for dark matter search experiments. 

Here, we should note that the Q-ball decay would produce too large amount of dark 
matter density if there were no pair annihilation of the LSP's. This can be easily seen 
by integrating the Eq. ( |T0"D with (av) = 0: 



Xlnoann \ 3 J J I ) B 



' ^fx^)(^)(°4)\ (12) 



100 GeV/ V h 



where m n ~ 1 GeV is the nucleon mass, and we have used the bound on the present 
baryon density VLsh 2 ^ 0.004 fll8|. Therefore, in the case of the bino-like LSP, the Q-ball 



formation is a serious obstacle for the AD baryogenesis. (Detailed discussion on this 
problem and possible solutions are given in Ref. ||11||.) 



As we have seen, Higgsino-like and wino-like LSP are promising candidates for cold 
dark matter if the AD baryogenesis is responsible for generating the observed baryon 
asymmetry in the present universe. Encouragingly enough, if this is the case, the direct 
detecting possibility for these dark matter is enormously enhanced compared with the 
case of bino-like dark matter |19|, |13| . 

The relevant quantity for direct search experiments is the elastic neutralino-nucleon 



6 In Fig.pl the annihilation cross section of the Higgsino is dominated by the decay mode into ti 
because of the light stop when m x > mt- This is why the resultant relic abundance of the LSP seems 
almost constant contrary to the naive expectation. 
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scattering rate [pO]: 



ap haio v F 

R= Hx * \ 13 
m x M N v 1 

where p^ a/o and v x is the mass density and the average speed of the neutralinos in the 
galactic halo, respectively. is the mass of the target nucleus, and F% is the nuclear 
form factor. By using typical values, this scattering rate is written as 

« = 4-lix 10" GeV ( 5", \ ( ;\ ) (lOp-) . (14) 
m x M N V0.3 GeV/cm 3 / V 320km /secy V k g- da y/ 



The list of relevant coupling constants are given in Ref. |p0| . For the numerical calcu- 
lations in this work, we have neglected the squark and Z-boson exchange contributions, 



since they are subdominant components in most of the parameter space fll3fl ) and we 
have taken F% — 1 for simplicity. 

In Fig. |] and |5], we show the scattering rates for the Higgsino-like and wino-like 
dark matter in 76 Ge detector, respectively. In these calculations, we have taken the 
same parameter sets as in Fig. [I] and except tan (3. From these figures, we see that the 
detection rates for the Higgsino-like and wino-like dark matter are R > 0.01 events/kg • 
day in very wide range of parameter space, and they even reach R > 0.1 events/kg • day 
in the large tan/? region. We should stress that such a parameter region where R > 
0.01 (events/kg ■ day) is within the reach of the on-going cold dark matter searches ||21|| . 

For comparison, we also show the scattering rate in 76 Ge detector for the case of 
bino-like dark matter in Fig. || Here, we have taken Mi = (l/3)/z, M<i = (2/3)/z, and 
other parameters are the same as in Fig. f| and |5|. As we noted, the scattering rate for 
the bino-like dark matter is much smaller than Higgsino-like and wino-like dark matter, 
and the direct detection in the dark matter search experiments is much more difficult. 



To summarize, we pointed out in this letter that the Higgsino-like and wino-like 
LSP's can be excellent dark matter candidates if the Affleck-Dine baryogenesis is re- 
sponsible for the generation of the observed baryon asymmetry in the present universe. 
Actually, we showed that the relic abundances of these LSP's can naturally explain the 
observed dark matter density with natural Q-ball decay temperatures, even in the rel- 
atively light neutralino mass region, which is much advantageous for the direct dark 
matter searches fZI . 



The novel thermal history of the universe proposed in this letter may have impor- 
tant implications on general SUSY breaking models, which include mSUGRA, no-scale 
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type models with non-universal gaugino masses p2| |, anomaly mediated SUSY breaking 



model | 23| , and so on. Detailed analysis on specific models will be given elsewhere |24 

In the context of the anomaly-mediated SUSY breaking (AMSB) [ ^3|| , non-thermal 
production of the wino dark matter from the late time decay of moduli and gravitino was 
investigated in Ref. [113. In the case of moduli decay, however, there exists large entropy 
production which substantially dilutes the primordial baryon asymmetry. The authors 
suggested that the AD baryogenesis can produce enough baryon asymmetry even in this 
case. Although the AD mechanism in AMSB scenario is generally difficult, there is an 
attractive AD scenario 1 1 L| which naturally works even in AMSB. However, if there is 
a large entropy production from the moduli decay, it is highly difficult to generate the 
required baryon asymmetry. Even if this is possible, the decay of the resultant large Q- 
ball probably plays a comparable role with the moduli decay as the non-thermal source 
of the LSP's. 
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Figure 1: The contour plots of the relic abundance Q x of the Higgsino-like neutralino 
LSP in the m x -Td plane. We have taken h = 0.7, Mi = (3/2)//, M 2 = 3//, tan/3 = 5, 
m^o = 300 GeV, a = 0, and mo = 1 TeV. The three shaded regions correspond to the 
range of Q x < 0.1, 0.1 < Q x < 1, 1 < Q x , from the top to the bottom, respectively. 
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Figure 2: The contour plots of the relic abundance Q x of the Higgsino-like neutralino 
LSP in the m x -Td plane. The parameters are the same as Fig. [I] except mo = 330 GeV. 
The three shaded regions correspond to the range of Q x < 0.1, 0.1 < Q x < 1, 1 < Q x , 
from the top to the bottom, respectively. 
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Figure 3: The contour plots of the relic abundance Q x of the wino-like neutralino LSP 
in the m x -T d plane. We took M\ = (3/2)// and M 2 = (1/2)//. Other parameters are 
the same as Fig. [I]. The three shaded regions correspond to the range of Q x < 0.1, 
0.1 < fL. < 1, 1 < Q x , from the top to the bottom, respectively. 
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Figure 4: Contours of the detection rate for the Higgsino-like dark matter in 76 Ge de- 
tector. The four shaded regions correspond to the ranges of the detection rate R > 0.1, 
0.1 > R > 0.03, 0.03 > R > 0.01, 0.01 > R (events/kg ■ day) from left to right, 
respectively. The parameters used in this calculation are the same as in Fig. [I] except 
tan (3. 
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Figure 5: Contours of the detection rate for the wino-like dark matter in 76 Ge detector. 
The four shaded regions correspond to the ranges of the detection rate R > 0.1, 0.1 > 
R > 0.03, 0.03 > R > 0.01, 0.01 > R (events/kg • day) from left to right, respectively. 
The parameters used in this calculation are the same as in Fig. |3] except tan (3. 
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Figure 6: Contours for the detection rate for the bino-like dark matter in w Ge detector. 
The four shaded regions correspond to the ranges of the detection rate 0.03 > R > 0.01, 
0.01 > R > 0.003, 0.003 > R > 0.001, 0.001 > R (events/kg • day) from left to right, 
respectively. Here we have taken Mi = (l/3)/x and M 2 — (2/3)/x. Other parameters are 
the same as in Fig. |] and |5]. 
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